Unusual localisation effects in quantum percolation 
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We present a detailed study of the quantum site percolation problem on simple cubic lattices, 
thereby focussing on the statistics of the local density of states and the spatial structure of the 
single particle wavefunctions. Using the Kernel Polynomial Method we refine previous studies 
of the metal-insulator transition and demonstrate the non-monotonic energy dependence of the 
quantum percolation threshold. Remarkably, the data indicates a "fragmentation" of the spectrum 
into extended and localised states. In addition, the observation of a chequerboard-like structure of 
the wavefunctions at the band centre can be interpreted as anomalous localisation. 
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Disordered structures attracted continuing interest 
over the last decades, and besides the Anderson local- 
isation problem^ quantum percolation^''^ is one of the 
classical subjects of this field. Current applications con- 
cern e.g. transport properties of doped semiconductors^ 
and granular metals^, metal-insulator transition in t"wo- 
dimensional n-GaAs heterostructures^, 'wave propagation 
through binary inhomogeneous media^, superconductor- 
insulator and (integer) quantum Hall transitions^iS, or 
the dynamics of atomic Fermi-Bose mixturesiS. Another 
important example is the metal-insulator transition in 
perovskite manganite films and the related colossal mag- 
netoresistance effect, 'which in the meantime are believed 
to be inherently percolativepii 

In disordered solids the percolation problem is char- 
acterised by the interplay of pure classical and quantum 
effects. Besides the question of finding a percolating path 
of "accessible" sites through a given lattice the quantum 
nature of the electrons imposes further restrictions on 
the existence of extended states and, consequently, of a 
finite DC-conductivity. As a particularly simple model 
describing this situation we consider a tight-binding one- 
electron Hamiltonian, 



N 



H^J2 ^»^1^» -^J2(4cj + H.c. 



(1) 



on a simple cubic lattice 'with N ~ L^ sites and random 
on-site energies e^ dra"wn from the bimodal distribution 



p{ei) ^ p6{ei - e^) + (1 -p)(5(e,; - es) . 



(2) 



The t'wo energies ca and cb could, for instance, represent 
the potential landscape of a binary alloy ApBi_p, -where 
each site is occupied by an A or B atom -with probability p 
or 1— p, respectively. In the limit A = (e_B — sa) — ^ oo the 
'wavefunction of the A sub-band vanishes identically on 
the B-sites, making them completely inaccessible for the 
quantum particles. We then arrive at a situation "where 
non-interacting electrons move on a random ensemble of 
lattice points, "which, depending onp, may span the entire 



lattice or not. The corresponding Hamiltonian reads 

H = ~tJ2(cUj+'a-c.), (3) 

{ij)£A 

"where the summation extends over nearest-neighbour A- 
sites only and, "without loss of generality, eA is chosen to 
be zero. 

Within the classical percolation scenario the percola- 
tion threshold Pc is defined by the occurrence of an infi- 
nite cluster Aoo of adjacent A sites. For the simple cubic 
lattice this site-percolation threshold is Pc = 0.3117ii^ In 
the quantum case, the multiple scattering of the particles 
at the irregular boundaries of the cluster can suppress the 
"wavefunction, in particular, "within narro"w channels or 
close to dead ends of the cluster. Hence, this type of dis- 
order can lead to absence of diffusion due to localisation, 
even if there is a classical percolating path through the 
crystal. On the other hand, for finite A the tunnelling be- 
t"ween A and B sites may cause a finite DC-conductivity 
although the A sites are not percolating. Naturally, the 
question arises "whether the quantum percolation thresh- 
old pg, given by the probability above "which an extended 
"wavefunction exists "within the A sub-band, is larger or 
smaller than pc- Previous resultsi^ for finite values of A 
indicate that the tunnelling effect has a marginal influ- 
ence on the percolation threshold as soon as A ::^ AtD, 
"where D denotes the spatial dimension of the hypercubic 
lattice. 

In the theoretical investigation of disordered systems 
it turned out that distribution functions for the ran- 
dom quantities take the centre stageii^ The distribution 
f{pi{E)) of the local density of states (LDOS), 
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p,(i?) = ^|^„(rOP<5(ii;-S„), 



(4) 



n=l 



is particularly suited because Pi{E) measures the local 
amplitude of the "wavefunction at site r^. It therefore 
contains direct information about the localisation prop- 
erties. In contrast to the (arithmetically averaged) mean 
DOS, Pmc{E) = {pi{E)), the LDOS becomes critical at 



the localisation transitiori j^^i^^ The probability density 
f{pi{E)) was found to have essentially different proper- 
ties for extended and localised phases 41 For an extended 
state at energy E the amplitude of the wavefunctions 
is more or less uniform. Accordingly f{pi(E)) is sharply 
peaked and symmetric about pmc{E). On the other hand, 
if states become localised, the wavefunction has consid- 
erable weight only on a few sites. In this case the LDOS 
strongly fluctuates throughout the lattice and the corre- 
sponding LDOS distribution is very asymmetric and has 
a long tail. Above the localisation transition the distri- 
bution of the LDOS is singular, i.e., mainly concentrated 
at Pi = 0. Nevertheless the rare but large LDOS-values 
dominate the mean DOS pmc{E), which therefore can- 
not be taken as a good approximation of the most prob- 
able value of the LDOS. Such systems are referred to 
as "non-self-averaging" . Of course, for practical calcula- 
tions the recording of entire distributions is a bit incon- 
venient. Instead the mean DOS Paic{E) together with 
the (geometrically averaged) so-called "typical" DOS, 
Pty{E) = exp(lnpi(£')), is frequently used to monitor the 
transition from extended to localised states. The typical 
DOS puts sufficient weight on small values of pi and a 
comparison to PmciE) therefore allows to detect the lo- 
calisation transition. This has been shown for the pure 
Anderson modeUSiiSiiSi and for even more complex situa- 
tions, where the effects of correlated disorder^", electron- 
electron interaction2L22, or electron-phonon coupling22iS4 
were taken into account. 

In this paper we employ the typical-DOS concept to 
analyse the nature of the eigenstates (extended or lo- 
calised) of the Hamiltonians ^ and (|3J) . Using the Ker- 
nel Polynomial Method-'^^'^^, an efficient high-resolution 
Chebyshev expansion technique, in a first step, we cal- 
culate the LDOS for a large number of samples, Kr, and 
sites. Kg. The mean DOS is then simply given by 



Pmc{E) 
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(5) 



'fc=l i=l 



whereas the typical DOS is obtained from the geometric 
average 



1 



K,. K, 



pty{E) = exp Tr^EEl^^('^^(^)) 



' fe=i 1=1 



(6) 



We classify a state at energy E with pmc(-E) 7^ as lo- 
calised if Pty{E) = and as extended if Pty{E) ^ 0. 

Before discussing possible localisation phenomena let 
us investigate the behaviour of the mean DOS for the 
quantum percolation models (Q and l|21l. Figure ^ shows 
that as long as ea and e^ do not differ too much there 
exists an asymmetric (if p ^ 0.5) but still connected elec- 
tronic band-i^ At about A ~ AtD this band separates 
into two sub-bands centred at e^ and es, respectively. 
The most prominent feature in the split-band regime is 
the series of spikes at discrete energies within the band. 
As an obvious guess, we might attribute these spikes to 
eigenstates on islands of ^ or i? sites being isolated from 
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FIG. 1: (Colour online) Mean (solid line) and typical (dashed 
line) DOS of the Hamiltonian (0 on a 50"^ lattice with periodic 
boundary conditions (PBC). Results are obtained using M — 
32768 Chebyshev moments and KaXKr = 32x32 realisations. 
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FIG. 2: (Colour online) Mean (solid line) and typical (dashed 
line) DOS for the quantum percolation model in the limit 
A — + 00. While in the upper left panel all A-sites are taken 
into account, the other three panels show data for the re- 
stricted model Q on the spanning cluster Aao only (note 
that pty is smaller in the former case because there are more 
sites with vanishing amplitude of the wavefunction) . System 
sizes were adapted to ensure that Aoo always contains the 
same number of sites, i.e., A^ — 57"^ for p = 0.405, 46^ for 
p = 0.70, and 42^ for p = 0.92. Again we used M = 32768 
and Ks X Kr = 32x 32. 



the main cluster,^'^^ It turns out, however, that some of 
the spikes persist, even if we neglect all finite clusters 
and restrict the calculation to the spanning cluster of A 
sites, Aoo ■ This is illustrated in the upper panels of Fig.|21 
where we compare the DOS of the models (^ [at A ^ cx)] 
and ^. Increasing the concentration of accessible sites 
the mean DOS of the spanning cluster is evocative of the 
DOS of the simple cubic lattice, but even at large values 
of p a sharp peak structure remains a,t E = (cf . Fig. |21 
lower panels). 

To elucidate this effect, which partially is not ac- 
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FIG. 3; (Colour online) Left: Mean (solid line) and typical 
(dashed line) DOS for the model with p = 0.33 on a 100^ 
lattice (PBC, M — 32768). Data obtained from an average 
over 100 random initialisations of sites for K^ = 100 realisa- 
tions of disorder. Note that pty < 10~^ holds in the whole 
band. Right: Some cluster configurations related to the spe- 
cial energies at which the peaks in pme occur. 



counted for in the literature2iii2LS&, in more detail, 
in Fig. 121 we fixed p at 0.33, shortly above the clas- 
sical percolation threshold, and increased the ensem- 
ble size. Besides the most dominant peaks at E/t ~ 
0, ±1, ±\/2, in addition we can resolve distinct spikes at 
E/t = i (±1 ± 75) , ±\/3,±\/2±y2, . . . These special 
energies coincide with the eigenvalues of the tight-binding 
model on small clusters of the geometries shown in the 
right part of Fig. O In accordance with Refs. '2' and 1291 
we can thus argue that the wavefunctions, which cor- 
respond to these special energies, are localised on some 
"dead ends" of the spanning cluster. 

The assumption that the distinct peaks correspond to 
localised wavefunctions is corroborated by the fact that 
the typical DOS vanishes or, at least, shows a dip at these 
energies. Occurring also for finite A (Fig.^l, this effect 
becomes more pronounced as A — > cx) and in the vicinity 
of the classical percolation threshold Pc- From the study 
of the Anderson modeli we know that localisation leads 
at first to a narrowing of the energy window containing 
extended states. The corresponding mobility edges have 
been mapped out with high precision^Siffi For the perco- 
lation problem, in contrast, with decreasing p the typical 
DOS indicates both localisation from the band edges and 
localisation at particular energies within the band. Since 
finite cluster wavefunctions like those shown in Fig.|31can 
be constructed for numerous other, less probable geome- 
tries^, Chayes et al.— argued that an infinite discrete 
series of such spikes might exist within the spectrum. 
The picture of localisation in the quantum percolation 
model is then quite remarkable. If we generalise our nu- 
merical data for the peaks at £^ = and E/t = ±1, 
it seems as if there is an infinite discrete set of energies 
with localised wavefunctions, which is dense within the 
entire spectrum. In between there are continua of de- 
localised states, but to avoid mixing, their density goes 
to zero close to the localised states. Facilitated by the 
large special weight of the peak (up to 10% close to Pc) 
this is clearly observed at E = 0, and we suspect similar 
behaviour at E/t = ±1. For the other discrete spikes the 
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FIG. 4: (Colour online) Amplitudes of the wavefunction, 
|V'n(ri)|, of the quantum percolation model Q on Aoo of a 14"^ 
lattice with occupation probability p = 0.45. Exact diagonal- 
isation results are presented at three characteristic energies 
E/t — 0, 1, and 0.66 corresponding to anomalously localised, 
localised, and delocalised states, respectively. 



resolution of our numerical data is still too poor and the 
system size might be even too small to draw a definite 
conclusion. 

In order to understand the internal structure of the 
extended and localised states we calculated the ampli- 
tudes of the wavefunction at specific energies for a ran- 
dom sample of the quantum percolation model ^ re- 
stricted to Aoo- Figure 0] visualises the spatial variation 
of |f/'„(ri)| on a 14'^ lattice with PBC and an occupa- 
tion probability p — 0.45 well above the classical percola- 
tion threshold. The figure clearly indicates that the state 
with E/t = 0.66 is extended, i.e. the spanning cluster is 
quantum mechanically "transparent" . On the contrary, 
a.t E — t, the wavefunction is completely localised on a 
finite region of the spanning cluster. Here the scatter- 
ing of the particle at the random surface of the spanning 
cluster results in states, where the wavefunction vanishes 
identically except for some finite domains on loose ends 
(like those shown in Fig. ^, where it takes the values 
(±1,=f1), (±1,=f1,0,t1,±1), • ■ • Note that these regions 
are part of the spanning cluster, connected to the other 
sites by a site with wavefunction amplitude zero. A par- 
ticularly interesting behaviour is observed ai E = 0. The 
eigenstate E — is highly degenerate and we can form 
wavefunctions that span the entire lattice in a chequer- 
board structure with zero and non-zero amplitudes (see 
Fig. 2J| . Although these states are extended in the sense 
that they are not confined to some region of the cluster, 
they are localised in the sense that they do not contribute 
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FIG. 5: (Colour online) Amplitudes |-!/),i(ri)| of a localised 
state on a 21^ lattice with p — 0.33 obtained by exact diago- 
nalisation of the model 0. E/t = 3.456 was chosen in order 
to avoid any of the special cluster configurations discussed 
above. 



to the DC-conductivity. This is caused by the alternating 
structure which suppresses the nearest-neighbour hop- 
ping, and in spite of the high degeneracy, the current 
matrix element between different E = states is zero. 
Hence, having properties of both classes of states these 
states are called anomalously localised4i2^ The chequer- 
board structure is also observed for the hypercubic lattice 
in 2D but with reduced spectral weight, compared to the 
3D case. Another indication for the robustness of this 
feature is its persistence for mismatching boundary con- 
ditions, e.g., periodic (antiperiodic) boundary conditions 
for odd (even) values of the linear extension L. In these 
cases the chequerboard is matched to itself by a manifold 
of sites with vanishing amplitude. 

In the past most of the methods used in numerical 
studies of Anderson localisation have also been applied 
to the percolation models ^ or ||2Jl in order to deter- 
mine the quantum percolation threshold pq, defined as 
the probability p below which all states are localised (see, 
e.g., Refs. lalSa and references therein). However, so far 
the results for pq are far less precise than, e.g., the values 
of the critical disorder reported for the Anderson model. 
For the simple cubic lattice numerical estimates of quan- 
tum site-percolation thresholds range from 0.4 to 0.5. In 
Figs. I1I3I we presented data for pty which shows that 
Pq > Pc- In fact, within numerical accuracy, we found 
Pty — ioT p = 0.33 > Pc- Figure [S] displays the ampli- 
tude of a typical state obtained by exact diagonalisation 
for a fixed realisation of disorder. Bearing in mind that 
we have used PBC, the support of the wavefunction ap- 
pears to be a finite (connected) region of the spanning 
cluster Aao, i-C, the state is clearly localised. 

To get a more detailed picture we calculated the nor- 
malised typical DOS, R{p, E) — pty/Pmc, in the whole 
concentration-energy plane. Figure El presents such kind 
of phase diagram of the quantum percolation model ^ . 
The data supports a finite quantum percolation thresh- 
old Pq ^ 0-4 > Pc (cf. also Refs. .13-33-34.35.) . but as the 
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FIG. 6: (Colour online) Normalised typical DOS R — pty/pmc 
in the concentration-energy plane for the model lO on a 50'^ 
lattice for p > 0.5 and a 100^ lattice for p < 0.5. We used 
M — 16384 Chebyshev moments and averaged over Ks x Kr = 
32 X 32 realisations. 



discussion above indicated, ioi E ~ and E = ±t the 
critical value Pq{E) is 1, and the same may hold for the 
set of other "special" energies. The transition line be- 
tween localised and delocalised states, Pq{E), might thus 
be a rather irregular (fractal?) function. On the basis 
of our numerical treatment, however, we are not in the 
position to answer this question with full rigour. 
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FIG. 7: (Colour online) Characteristic probability distribu- 
tions of the LDOS. 

Finally let us come back to the characterisation of ex- 
tended and localised states in terms of distribution func- 
tions. Figure [7| displays the probability distribution. 



F[P^{E)] = 



P.{E) 



f{p',{E))dp[{E), 



(7) 



for three typical energies E/t = 3.5, 1.6, and cor- 



responding to localised, extended, and anomalously lo- 
calised states at p = 0.45, respectively. The differences 
in F[pi{E)\ are significant. The slow increase oi F[pi{E)\ 
observed for localised states corresponds to an extremely 
broad LDOS-distribution, with a very small most proba- 
ble (or typical) value oi pi{E). This is in agreement with 
the findings for the Anderson model. ^^-^^ Accordingly the 
jump-like increase found for extended states is related to 
an extremely narrow distribution of the LDOS centred 
around the mean DOS, where Pnic coincides with the 
most probable value. At ii^ = 0, the probability distribu- 
tion exhibits two steps, leading to a bimodal distribution 
density. Here the first (second) maximum is related to 
sites with a small (large) amplitude of the wavefunction 
- a feature that substantiates the chequerboard structure 
discussed above. 

To summarise, we have demonstrated the value and 
power of the probability distribution approach to quan- 
tum percolation. As for standard Anderson localisation 



the typical density of states can serve as a kind of order 
parameter differentiating between extended and localised 
states. Our numerical data corroborates previous results 
in favour of a quantum percolation threshold Pq > Pc 
and a fragmentation of the spectrum into extended and 
localised states. The latter refers to a discrete but dense 
set of localised states separated by continua of delocalised 
states. Accordingly the function Pq{E) is rather irregu- 
lar. At the band centre, so-called anomalous localisation 
is observed, which manifests itself in a chequerboard-like 
structure of the wavefunction. Even though the Kernel 
Polynomial Method allows for the study of very large 
clusters with high energy resolution, the quantum perco- 
lation problem certainly deserves further investigation. 
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